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The  important  collisional  energy  transfer  ratio,  j  =  <flE>a j j /<AE>j  was 

a 

examined  by  using  a  fitted  classical  approximation  for  the  density  of  internal 
energy  eigenstates.  The  exponential  form  of  the  collisional  energy  transfer 
probability  function  was  applied  to  four  model  unimolecular  reaction  systems. 
Two  parameters,  the  inversion  temperature  Tj  and  the  effective  temperature  T@ 
(defined  previously  in  I),  were,  employed  to  develop  a  form  for  the  parametric 

• '  ^ u  +'  -in¬ 

dependence  of  7  on  <AE>.j,  energy  level  E  and  temperature  T;  these  quantities 

•  y.-,  b  -  . 

are  related  to  a  reduced  average  energy  transfer  quantity  E^-.  Comparison  was 

A 

made  with  previous  literature  expressions  for  the  inter-relationship  between 

£  /  /  /•  O  i  *  -  -r  r  '  — 

<AE>ajj  and  <£E>^  for  both  exponential  and  stepladder  models  of  the 


SL 


Introducti on 


In  a  recent  paper*  (I),  we  discussed  the  parameterization  of  unimolecular 

reaction  rates  involving  weak  collisions  of  reactant  molecules  with  a  bath 

species.  The  collisional  vibrational  energy  transfer  efficiency  (3 can  be 

related  to  either  of  two  average  energy  quantities  conventionally  used, 

namely,  <AE>^  and  <AE>a|j.  The  former  is  the  collisional  internal  energy 

change  averaged  over  the  down  transitions,  while  the  latter  is  the  step  size 

averaged  over  both  down  and  up  transitions.  Both  quantities  have  been  widely 

2  3 

used  in  the  literature,  most  recently  by  Barker  and  Golden  (BG) .  Oref  has 

treated  these  quantities  for  the  case,  especially,  of  the  Boltzmann  model  of 

collisional  transition  probabilities. 

4 

Gilbert  has  recently  pointed  out  that  <AE>^  is  the  more  appropriate 

quantity  to  use  in  parameterizing  indeed,  this  is  the  quantity  which  is 

implicit  in  the  original  L i ndemann-Hi nshel wood  condition  on  collisional 

efficiency  in  their  formulation  of  unimolecular  reactions.  Although  the 

<AE>a||  quantity  is  directly  accessible  from  even  limited  experiment,  it  has 

5 

limited  physical  content.  For  good,  comprehensive  data,  however,  one  can 
also  extract  (AE)^  which  is  required  for  the  specification  of  the  correct  form 

g 

of  the  collisional  transition  probability  function.  Theoretical  progress,  in 
effect,  is  closely  coupled  to  the  further  elucidation  of  this  form. 

It  may  be  pointed  out  that  the  subject  of  weak  gas  collisions  and  their 
role  in  unimolecular  reactions  is  not  of  relatively  recent  origin  as  is  sometime 
supposed.  These  questions 

were  entertained  first  in  the  late  twenties,  while  in  the  thirties  O.K.  Rice 
and  D.V.  Sickman  made  quantitative  weak  collider  calculations.  In  the  early 
fifties,  following  the  hiatus  of  World  War  II,  H.S.  Johnston  and  coworkers 
re i nv i gorated  the  whole  question  of  the  pressure,  energy,  and  species  depen¬ 
dence  of  the  .trength  of  collisions.  A  detailed  review  may  be  found  in  ref  5. 


In  I,  we  def  ined  the  quantity  7  as  the  ratio  <AE>a  j  j /<AE>cj .  The 

magnitude  of  7  depends  strongly  on  the  form  of  the  coilisional  transition 

probability  matrix  P  that  is  appropriate  (or  selected).  We  may  remind  the 

reader  that  the  step  I  adder  (SL)  and  the  exponential  (EXP)  models  for  the 

transition  probabilities  are  the  two  types  of  models  which  have  been  used  most 

prominently;  the  SL  model  may  be  considered  as  an  approximation  for  the  more 

physically  apt  gaussian  (G)  distribution.  [More  recently,  a  reverse 

2  4 

exponential  model  has  also  been  used  '  and  has  properties  similar  to  that  of 

the  SL  and  G  models,  although  any  deeper  physical  connotation  of  this  form 

seems  more  obscure.]  The  SL  and  G  models  correspond  to  the  physical  case 

where  large  downjump  transitions  from  the  initial  energy  level  are  more 

probable  than  are  very  small  transitions.  While  the  EXP  model  corresponds  to 

the  reverse  situation,  where  very  small  transitions  have  higher  probability 

5  7 

than  larger  transitions.  It  has  been  shown  '  that  the  former  model  applies 
to  cases  where  the  behavior  corresponds  to  stronger,  more  efficient 
collisions,  and  that  the  EXP  model  is  more  correct  for  very  weak,  inefficient 
collisions.  As  a  rough  rule  of  thumb,  the  dividing  line  between  the  two 
models  of  behavior  corresponds  to  p  1  0.25. 

The  magnitudes  of  <AE>^  and  <AE>g | |  may  be  similar  or  quite  disparate, 
depending  upon  the  nature  and  conditions  of  a  particular  system,  so  that  a 
relationship  between  the  two  is  also  desirable  and  necessary  for  the  inter¬ 
pretation  of  literature  results  which  may  be  cast  in  either  form.  The  general 

and  correct  interconversion  between  these  quantities  is  not  facile.  A 

8  9 

relation  was  given  by  Troe  '  for  the  EXP  model  as, 

<  &E  >aii  =  <  A 2?  >2  /(<  A.E  >d  — F[-:RT )  , where  the  density  ratio  is 

Fe  =  \exp(E0/RT)/RT»(E0)\  f  p(E)exp(-E/RT)dE 

J  Ea 


was  emphasized  to  be  near  unity  ,  usually,  although  its  deviation  from 
9 

unity  was  also  shown.  In  fact,  departure  from  unity  is  important  under  many 
experimental  conditions. 

i  8  9 

BG  have  recently  given  a  critique  of  Troe  s  earlier  treatment0'  for  the 
EXP  model  and  offered  an  improved  correlation  expression  between  the  two 
energy  parameters  for  both  the  EXP  and  SL  cases.  They  used  the  Whitten- 
Rabinovitch  (WR)  approximation  to  give  an  expression  for  a  desired  vibrational 
eigenstate  density  ratio.  Notwithstanding  the  obvious  merit  of  the  BG 
equations,  it  is  very  desirable  to  test  other  approximations  that  interrelate 
and  rationalize  these  quantities.  An  iterative  analytical  procedure  was 
developed  in  I  for  interrelating  the  two  quantities  for  the  case  of  the  SL 
model.  In  this  treatment,  a  classical  expression  for  the  density  of 
vibrational  states  was  employed  in  which  the  ef f ecti ve  number  of  vibration 
modes  s' (E)  was  evaluated  from  the  correct  quantum  statistical  expression. 

Two  useful  temperature  parameters,  Tj  the  inversion  temperature,  and  Tg  the 
effective  temperature  were  also  introduced. 

In  the  work  described  below,  we  have  extended  the  comparison  given  in  I 
between  <AE>a||  and  <AE>^  to  the  EXP  model  case.  Comparison  of  both  models  is 
also  made  with  earlier  literature  expressions. 

[We  note,  in  response  to  a  referee,  that  although  the  discussion  here  is 
couched  in  terms  of  reactive  systems  and  illustrated  below  for  systems  at 
the  level  of  the  critical  threshold  energies,  the  conclusions  derived  are 
actually  more  general,  of  course,  and  apply  to  non-reactive  systems,  as  well, 
and  to  any  chosen  energy  level  (such  that  the  approximation  conditions 
incorporated  in  our  calculations  are  not  violated).] 


Ca leu  I ational 


Vibrational  densities.  Most  of  the  calculations  were  made  for  the  same 
prototype  reactions  as  in  I:  the  decomposition  of  nitryl  chloride,  and  the 
isomerizations  of  methyl  isocyanide,  cyclopropane  and  eye  I oheptatr i ene  (CHT) . 
These  reactants  reflect  differences  in  vibrational  frequency  patterns  and 
molecular  complexity,  excitation  levels  and  reaction  temperatures  (Table  I). 

The  energy  dependence  of  the  density  of  vibrational  states  (E  in  the 
classical  limit)  can  be  parameterized  either  by  adding  the  WR  energy  dependent 
correction  term  involving  the  zero  point  energy,  i.e.  [E  *  a(E)  E^]  ,  or  by 

taking  s  to  be  energy  dependent,  s' (E) ;  in  the  latter  case,  s' (E) 
monotonical ly  increases  with  energy  to  the  maximum  number  of  oscillators,  s 
(note  that  internal  rotations  may  be  easily  incorporated). 

Exact  calculations  of  reference  densities  as  bench  marks  for  all  of  the 
approximations  examined  here  were  made  using  densities  based  (arbitrarily)  on 
the  Haarhoff  algorithm.  No  differences  of  any  consequence  arise  in  any  of  the 
calculations,  figures,  tables  or  conclusions  presented  below  if,  instead,  WR 
densities  are  employed  as  reference  values. 

Collisional  transition  probabilities.  The  equations  representing  the 
down  probabilities  for  the  SL  and  EXP  models  are,  respectively, 
p$L  =  p(E, ,  Ej)  -  const.,  for  E}  -  E,  =  A£ 

-  0,  for  Ej  -  Et  A E 

and 

PdXP  =  P(E  Ej)  =  const,  exp  [  -  (E}  -  £,)/  <  A E  >d  ],  for  E,  <  Ej 

Up  transition  probabilities  p^  are  related  to  the  down  probabilities  p^ 
by  detailed  balance  and  completeness: 


p„(A£)/pd(A£)  =p(Ei,Ej)/p(EJ,Ei) 


=  p(Ei)/p{Ej)  exp  [  -  A E/RT]  ,  A E  =  \  Et  -  E,\  ;  (l) 

p(£i,  £j)  =  i 

all£. 

It  is  convenient  to  render  eq.  (1)  into  the  form, 

pu(A£)/pd(A£)  =  cxp(-A£/£Te)  ; 

which  requires  that  the  density  ratio  take  the  form  ,  p(Ej)/p(Ej)  =  exp(cAE) ; 

where  T  is  an  effective  temperature  given  by 
e 

Te  =  T/(l-  cRT)  (2) 

Table  1  illustrates  the  density  ratio  as  a  function  of  AE  for  various 
reactants,  at  Ej  =  Eq.  Of  prime  importance  is  the  near- I i near i ty  of  the  log  of 
the  density  ratio  with  increase  of  AE.  By  parameterizing  p(E)  as  ES  in 

I,  the  relation  £n(p(E+AE)/p(E))  =  (s'  (E)-l)£n(l«-AE/E)  was  used.  In  this 
paper,  we  wi I  I  call  the  fitted  classical  approximation,  designated  as  TR,  as 
the  s' (E)  value  which  makes 

ln[p{E+  <  AE  >d)/p(E)j  =  (s'(£)  -  l)  <  A E  >d  /E  (3) 

=  cjrAE  for  SL 


=  cTr  <  AE  >d  for  EXP,  with  <  AE  >d  /E  «  1. 
Alternatively,  BG  have  used  other  parameterization,  namely,  the  WR 
expression  whereby 


AE 


tn^A€)  =  (j  _  l  )fa  [£+«(£  ±  + 

p<f)  E  +  a(E)Ez  E  +  a(E)Ez  J 


(s  -  l)fn(l 


AE 


E  +  a(E)Et 


) 


(s-_l)A£  = 

E  +  a{E)Ez  w* 


(4) 


This  compact  formulation  has  been  made  to  depend  on  only  a  single  correction 
constant,  a(E). 


5 


Average  energy  quantities.  The  average  energy  quantities  may  be 

expressed  in  terms  of  the  transition  probabi I ities 

<  AJ E  >u  =  (Ej  —  Et)p(E},  E{)  j  P (Ej,E,) 

e3>ex  e,>e. 

<  A E  >a  =  —  (Ej  ~  Ei)p[Ej,Et)/  ^  p[Ej,E{) 

E,  <  E,  Ej  <  E, 

<  A£  >all  =  Pu  A E  >u  — Pd  A E  >d 

Where'  pu  =  P(EJ'E *)  and  Pd  -  X] 

These  expressions  take  an  obvious  simplified  form  for  the  SL  model. 


The  sign  convention  makes  <AE>a j  j  =  -<AE>d  when  p^  ■*  1,  and  <AE>a||  =  <AE>u 
when  p^  ■+  0.  Both  <AE>u  and  <AE>a | j  are  strongly  dependent  on  molecular 
complexity  and  temperature  for  a  given  <AE>d- 


The  temperature  quantities,  T^  and  Tj.  The  ratio  Pu/Pj  's  related  to  AE 
through  Tfi  which  is  independent  of  AE  for  c  constant.  The  condition  in  eq. 

(3)  is  that  <AE>d«  E;  in  eq.  4,  the  condition  is  <AE>d«  (E+a(E)Ez). 

For  molecules  of  given  energy,  e.g.  Eq,  an  inversion  temperature  Tj  was 
defined  in  I  as  the  temperature  at  which  <AE>g j |  =0  i .e.  <AE>d  =  <AE>u,  or, 

-  £  (Ey-£,)p(E1.E,)=  52  (E,-EME„E.)  (5a) 

Ej  <  Ei  Ej  >  £, 

For  this  constraint,  the  Boltzmann  factors  are  offset  by  the  density  ratios. 
With  a  stepladder  probability  model  at  energy  E  ,  this  constraint  was 
satisfied  (cf .  I)  when,  from  eq.  3, 

Ti  =  A E/tn{p{E0  +  A E)/p(E0))R  =  E0/(s’(E)  -  1)R . 

Consider  the  general  limiting  cases  for  Tg  (eq.  2): 

cRT  >>  1,  Te  =  —l/cR 


cRT  «  1,  Te  =  T 
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For  the  fitted  classical  approximation, 

.hen  ctrRT  =  1,  TjK  =  oo-,  so  r  =  7, ™  =£„/(*'(£)  -  1)»  l56) 

TR 

In  the  limit,  T  «  Tj  (i.e.  c-j-rRT  <<  1),  i.e.  at  low  temperatures  and/or  high 
TR 

E  ,  T  =  T  ;  in  this  limit,  p  <  p^,  the  maximum  in  the  Boltzmann  distribution 

TR 

is  below  Eo  and  down  transitions  dominate.  For  the  case,  T  >>  Tj  (CjrRT  » 

1) ,  i.e.  at  very  high  temperatures,  the  maximum  in  the  Boltzmann  distribution 

is  above  E  and  up  transitions  dominate. 
o  r 

We  may  now  derive  specific  expressions  from  these  general  results  for  the 
EXP  model  and  recall  expressions  from  I  for  the  SL  case.  The  expressions  are 
for  the  energy  origin  E-  =  Eo>  but  are  valid  for  any  E  so  long  as  <AE>^  <<  E 
for  the  EXP  model . 

Stepladdsr  model .  For  the  SL  model,  we  showed  in  I  that 


Pm!  PA  =  exp(-A  E/RTe) 


*7SL  =  -  (1-exp  '  -  A E'RT'j. 


1  -  exp(-A£ ;RT')\ 


—  —  tanh(A£’/2I2T'e)  (G) 

Thus  <AE>j  and  <AE>a j |  are  universally  related  via  Tg.  The  validity  of  the 
approximations  in  using  T^,  i.e.  in  using  a  constant  Tj,  is  illustrated  in 
Table  2  where  values  of  Tj  are  presented;  Tj  is  nearly  independent  of  step 
size:  there  is  a  slight  increase  in  Tj  with  an  increase  in  step  size  since 
the  increase  in  density  ratio  with  energy  does  not  completely  compensate  for 
the  decrease  in  the  Boltzmann  factor.  It  is  evident  that  a  complex  molecule 
(large  s)  with  a  large  Eq  will  have  a  comparable  Tj  to  a  smaller  molecule  with 
a  small  EqI  i.e.  (s  -1)/Eq  constant;  e.g.  C^Hg  and  CH^NC 
(Table  2) . 
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Exponential  model.  For  the  EXP  model 

<  A E  >aii=  }  " ( F  j  —  Ei)p{Ej  '  F\ ) 

E, 


Ze,<e,{E>-  -£,)/<  A£>d[  +  DE|>E|(£,  -  E.)cxpj-(£,  -  £,)/<  A£  >d]  *  fo,)exp[ -(£,  -  £,)//?T]/^ 


T.e,<e,  exp(-(£,  -  £>)/  <  A E  >d|  +  Ee,>e.  pxP!  ~(Ej  ~  £.)/  <  >dj  •  ^expj-(£,  -  E,)/RT |/p£,  ) 

Eae>0A£:1~1  -rexp(-A£//?rt)lexp(-A£/  <  A£  >d) 

1  +  Eae>oU  +  exp(-A£/i?re)|exp(-A£/  <  A£  >d) 

In  the  evaluation  of  eq.  7,  the  grain  size  G  (which  corresponds  to  the  minimum 
step  size)  was  taken  sufficiently  small,  and  the  number  of  steps  included  in 
the  EXP  model  sufficiently  large  so  that  convergence  was  obtained.  In  these 
calculations,  for  the  grain  size,  G  =  <AE>^/256,  and  an  energy  limit  of 
16<AE>j  (i.e.  4096  steps),  absolute  convergence  was  found  to  within  0.2%  in 
the  worst  case  and  was  <  0.1%  in  most  cases;  the  error  for  the  ratios  of  sums 

is  even  less.  In  the  limit,  AE  - ->  0,  the  sums  in  eq.  7  can  be  replaced  with 

an  integral  as  done  by  Troe  and  by  BG.  Integration  gives 


and 


<  A  £  >  ail  =  -  <  AE  >3  /(<  AE  >d  +RTe) 


=  -  <  a E>d  /(<  A E  >d  +RT™) 


where  TjR  =  17(1  -  (s'(£)  -  l)RTjE0),  since  TjR  =  T /(l  -  T jTjR). 


(8a) 

(86) 


8,9 


Earlier,  Troe  '  has  used, 

<  A E  > ail  =  -  <  A E  >3  /(<  A E  >d  +FeRT) 


so, 


7r*P  -  -  <  A E  :  (i  /(<  A E  >d  -FeRT) 

where  was  given  as  positive  and  was  evaluated  with  use  of  the  WR 
approximation. 


(9) 


(7 


BC  gave  the  relation, 

.  A  _  1  1  (S  -  1)  I"1 

>all“  +  RT  -  £T^£)£j  ~<AE><* 

We  may  simplify  this  equation  by  writing 


where, 


=  -  <  A E>d  /(<  A E  >d  +  RTeBG) 
tbg  =  Tj  ^  _  (5  _  i)RTj(E  +  a(E)E~)J  =  T/(l  -  T/T/30) 


EXP  T  T 

These  7  expressions  (eqs.  8a,  9  and  10)  are  equivalent  when  F^  =  T^1  = 

TR  BG 

T  =  T°  ;  this  condition  is  realized  for  small  molecules  at  low  temperatures. 
©  © 

The  general  behavior  is  illustrated  in  Table  3.  For  a  large  molecule  and 

higher  temperatures,  the  aberration  of  Fc  =  T^/T  from  T  (exact)/T  is  striking. 

t  0  e 

For  an  EXP  model,  eq.  5a  takes  the  form 


A E exp (  —  AiT/  <  A E  >a)  -  AE cxp(-  AE  <  A E  >d ) 


AE>0 


AE>0 


(p(p+ae)/p(f)  exp(- A E/RT) 


LAI 

under  the  condition,  T  =  Tj  .  This  equation  may  be  solved  iteratively  for 

EXP  EXP 

Tj  .  A  listing  of  values  is  given  in  Table  2.  It  is  noted  that  Tj  are 

somewhat  larger  than  those  calculated  for  the  SL  model  and  increase  somewhat 

faster  with  <AE>^  than  for  the  SL  model.  In  effect  a  higher  temperature  is 

required  to  offset  the  decrease  in  pu  due  to  the  higher  weighting  of  the  small 

step  sizes  (head  of  the  distribution)  characteristic  of  the  EXP  model. 


Results  and  Discussion 

SL  and  EXP  model  calculations  were  performed  for  the  four  reactions  over 

the  temperature  range  250  K  to  4000  K  with  average  step  sizes  ranging  from  100 

to  1600  cm  V  For  comparative  purposes,  ethyl  bromide  and  methyl 

eye  I ohepatr i ene  calculations  were  also  performed  for  two  temperatures  and  two 

step  sizes  each.  Results  for  these  compounds  are  in  good  agreement  with  those 

4 

published  by  Gilbert  .  Our  ethyl  bromide  results  at  1000  K  are  in  agreement 
with  the  corrected  values  of  Gilbert  (private  communication). 
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Step  ladder  model.  Results  are  summarized  in  Fig  1  where  7  is  plotted 

1  TR 

as  a  function  of  the  reduced  parameter,  E  ,  defined  as  <AE>  ,/RT  By 

e  d  6 

calibration  of  s' (Eq+AE)  for  every  AE,  the  slight  scatter  shown  in  Fig  2  of  I 

Cl 

has  been  removed;  the  calculations  follow  eq.  6.  As  expected,  7  increases 

from  -1,  for  large  <AE>^  and/or  low  temperature,  to  0,  for  small  <AE>^  and/or 

high  temperature,  and  continues  to  increase  for  temperatures  greater  than  Tj 

as  <AE>^  increases.  At  temperatures  greater  than  Tj,  up  transitions  dominate 

SL 

and  <AE>a)|  approaches  <AE>u  in  magnitude,  i.e.  7  =1.  For  most  small 

SL 

reactants,  the  region  where  7  >  0  is  not  experimentally  accessible.  For 

CHT,  Tj  was  realized  experimentally  by  Troe,  et.  al.^ 

Cl 

When  E'  <  0.2,  7  is  linear  with  E' ;  the  slope  is  -0.5.  This 
©  © 

observation  is  verified  by  expanding  eq.  (6)  which  gives, 

^SL  =  -A£  2 RTe  =  -  E'J 2 


Comparison  of  this  expression  with  the  limiting  form  of  eq.  8a,  gives  the 
EXP  SL 

relation  <AE>~7  =  AE,  /  \/2  ,  for  a  given  (experimental)  determination  of 


Exponential  Model.  The  results  for  the  EXP  model  are  illustrated  in 

Fig  1.  The  scatter  is  greater  than  for  the  SL  model,  but  the  fit  is 

remarkably  good.  The  deviations  are  due  to  the  large  energy  range  that  is 

EXP 

spanned  by  the  energy  transitions.  The  behavior  for  7  <0  is  qualitatively 

»  EXP 

similar  to  the  SL  model;  major  differences  occur  at  low  E^;  7  is  linear 

EXP 

with  E'  with  a  slope  of  -1  and  not  -0.5,  as  for  SL.  For  larger  E' ,  7 

SL  EXP 

approaches  -1  more  slowly  than  does  7  For  E^  <  0,7  exceeds  unity.  Thi 

results  from  the  fact  that  the  Boltzmann  temperature  factor  is  less  important 

than  the  density  ratio  increase;  larger  steps  are  enhanced. 

SL  EXP 

A  comparison  of  7  and  7  is  illustrated  in  Fig  2.  For  sufficiently 
large  E1 ,  <AE>^V.  =  <AE>Pf?  =  -<AE>,.  As  E'  decreases,  <AE>  ,,  for  the  EXP 


model  does  not  increase  as  fast  as  the  SL  model:  as  E*  ■*  0.  the  reverse 

e  ' 

becomes  true,  hence  the  minimum  in  Fig  2.  This  behavior  can  be  understood  by 

EXP 

looking  at  the  limiting  forms  of  eq.  7  incorporated  in  7 

For  small  E^,  the  quantity  exp(-AE/RTe)  may  be  expanded  in  linear  form 
and  the  summation  is  over  a  distrib:  V>on  which  depends  on  the  second  moment  of 
AE,  i.e.  favors  the  tail  of  the  energy  step  distribution  corresponding  to 
larger  steps.  When  this  sum  is  compared  to  a  SL  model  a  factor  of  2  results 
(the  limit  of  eq.  6) . 

For  large  E1 ,  exp(-AE/RT  )  is  close  to  0,  and  the  resulting  sum  in  eq.  7 

6  C 

has  only  a  linear  term  in  AE,  so  that  the  head  of  the  energy  step  distribution 

EXP 

dominates;  a  smaller  average  step  size  results.  For  E*  >  3,  7  approaches 

SL 

-1  more  slowly  than  does  7  ,  with  increasing  E^. 

We  note  that  the  behavior  in  Fig  2  is  reminiscent  of  the  variation  of  the 

EXP  SL  >  * 

relative  magnitudes  of  PQ  and  ^  taken  as  a  function  of  E1 (=<AE>d/<E*>) 

(a  parameter  similar  in  nature  to  E^)  which  was  described  by  us  some  time 

ago.^  Plots  of  sc  vs  E1  given  there  for  both  models  do  not  coincide  at 

low  values  of  E1  but  cross  and  converge  as  E1  increases.  This  behavior  can 

and  has  been  used  to  differentiate  transition  probability  models  and  their 

domain  of  relevance.  In  ref  11,  the  plots  of  Bc  vs  E‘  separate  by  a  factor 

of  ■v  /2  at  small  values  of  E'.  This  is  a  positive  feature  that  assists  in 

the  experimental  determination  of  the  correct  form  of  £,  and  is  not  a 

g 

shortcoming  as  misstated  by  Troe. 

Comparison  with  previous  work.  The  present  results  may  be  compared  to 

Q 

the  7  functions  presented  by  Troe  and  BG.  As  pointed  out  by  BG,  Troe' s 
function  (eq.  9)  suffers  from  the  qualitative  defect  that  it  cannot  provide 
7  P  >  Both  7EXJr  (eq.  8a)  and  7EXPQ  (eq.  10)  allow  for  7EXP  >  0,  since 
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can  be  less  than  zero.  This  behavior  may  be  analyzed  in  more  detail  as 
follows.  Equations  8a,  9  and  10  can  be  recast  into  the  form 


=  -£'/(i  +  E’e) 

For  eq.  8a,  E'c  =  <  A E  >d  / RT ™ ,  and  TjR /T  =  1/(1  -  T jTjR)\ 

for  eq.  10,  E'e  =<  A E  >d  /RTBG,  and  TeBG ,'T  =  1/(1  -T/Tf0); 

and,  for  eq.  9,  E't  =<  A E  >d  jRTj ,  and  Tj  >T  =  FE 


Tp  pg 

For  sufficiently  low  temperature,  the  quantities  T  /T,  T  /T  and  Fp  are 

TR 

all  nearly  linear  with  temperature.  At  higher  temperatures,  only  Tg  /T  and 
RP 

T°  /T  display  similar  coefficients  of  higher  powers  of  T.  Also,  when 

T  TR  BP 

T  >  Tj,  T^/T  remains  positive  while  Tg  /T  and  T°  /T  become  negative  as 

required  physically.  A  tabulation  of  the  Te/T  quantities  is  given  in 

Table  3  for  the  conditions  of  the  present  calculations. 

The  relative  merits  of  the  various  Tg/T  quantities  for  parameterizing 

EXP  .  . 

7  can  be  found  by  defining 


A^exp 

-yEXP (£/)  _  '-)(EXP (exact)  1 

^EXP 

7EXP(exact) 

The  resulting  values  of  R^  for  the  various  approximations  are  displayed  in 
Fig.  3  for  two  systems  (CHgNC  and  CHT)  over  a  range  of  temperature  and  step 
sizes,  although  calculations  were  made  for  all  four  molecules.  To  avoid 
display  of  three  scales  for  the  abscissa,  the  comparison  of  the  approximations 
is  shown  for  the  same  values  of  E^(TR),  although  the  appropriate  value  of 
E^(BG)  and  E^(T)  were  used  in  calculating  and  7j^  from  eq.  11a.  The 

comparison  is  correct  since  points  related  vertically  in  the  figure  do 
correspond  to  the  same  values  of  T,  <AE>  and  given  molecule.  To  be  noted  is 


the  increasing  deviation  of  R  for  all  three  approximations  as  E'  decreases- 

/  e 

this  occurs  for  either  a  decrease  in  step  size  and/or  an  increase  in 

temperature.  For  a  given  E^,  increases  with  molecular  complexity 

(CHT  y  CH^NC) .  For  all  calculations,  the  R^  deviations  are  largest  for  Troe's 

approximation-  hence  eq.  9  with  F^  is  not  recommended.  TR  values  are  better 

than  are  BG,  and  this  improvement  is  only  slightly  reduced  if  s' (E)  were  to  be 

used  as  a  constant  average  value  independent  of  <AE>. 

For  the  case  of  the  EXP  model,  the  BG  calculations  require  the  fewest 

number  of  function  evaluations.  Because  of  the  large  experimental  error  in 

measured  quantities,  <AE>^  or  <AE>g |  j ,  we  advocate  that  the  BG  approximation 

EXP 

be  used  unless  greater  ca Icul ationa I  precision  is  desired,  in  which  case  7 
(exact)  may  be  computed  readily  (indeed,  all  the  calculations  for  this  paper 
were  performed  on  an  IBM-PC) 

An  interesting  comparison  may  be  added.  Both  Troe  and  BG  derived  eqs.  9 
and  10  by  assuming  a(E)  =  a(EQ);  this  approximation,  although  not  generally 
justified,  is  appropriate  for  small  molecules  with  E  ~  Eq.  Troe  then 
performed  an  analytical  integration  to  determine  Fg(  =  T^/T);  BG,  on  the  other 
hand,  before  integration  introduced  a  fortunate  second  approximation 
£.n(l  ♦  AE/  (Eq  ♦  a(EQ)E)  a  AE/(Eq  ♦  a(EQ)Ez)  with  result: 

/  E„  -r  a{E)Ez  AE  \  /  AE  \  AE 

in\E0^-  a{E0)Ez  ~  E0  -r  a(E0)Ez  )  ~£nl1+  E0-ra[E0)Ej  ~  E0-ra{E0)Ez' 

and  with  near-compensation  of  errors  introduced  by  the  first  approximation. 
(See  Appendix). 

SL 

A  comparison  of  7  calculations  analogous  to  the  above  between  TR  and  BG 

is  superfluous  for  <5-function  transitions.  In  effect,  our  fitted  classical 

approximation  was  made  exact  for  any  initial  energy  Ej . 

Interrelation  of  <AE>,  and  <AE>  ...  The  transformation  from  <AE>,  to 
- d — 1 - a  I  I  0 

<AE>a | |  is  straight  forward:  i)  Tj  is  first  computed  by  eq.  5b;  ii)  Tg  is 

evaluated  for  the  specific  temperature  (eq.  lib);  iii)  E'  is  calculated  from 

6 
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<AE>d/RTe  and  used  in  eq.  6  or  eq.  11a  to  solve  for  7SL  or  7^  .  Finally, 

<AE>aj|  =  7  •  <AE>d  for  either  case. 

The  reverse  process,  to  derive  <AE>,  from  <AE>  ..  for  the  EXP  case, 

d  'all  ' 

utilizes  the  expression  (from  eq.  11a) 


<  A E  >j  = 


-  <  A E  >a|| 


1  +  v  1  -  4i?re/  <  A E  >a,i 


for  the  EXP  case.  As  described  in  I,  an  iteration  technique  is  used  to 
determine  <AE>d  for  the  SL  model. 

CVD 

At  the  request  of  a  referee  we  are  Including  Figure  4  where  y  (exact)  Is 

plotted  as  a  function  of  temperature  for  the  methyl  Isocyanide  and 

cyloheptatrlene  systems.  Step  sizes  of  100  and  1600  cm”1  are  Illustrated. 

EXP 

Values  of  y  (exact)  for  Intermediate  step  sizes  and/or  other  reactants  can  be 
found  by  Interpolation.  These  family  of  curves  Illustrate  the  advantage  of 
Incorporating  molecular  complexity,  step  size  and  temperature  In  the  reduced 

parameter  E1  so  that  a  quasl-unlversal  curve  results. 

Cone  I  us  ion 


A  critique  has  been  given  of  the  various  parameters  that  enter  into  the 

formulation  of  collisional  efficiency  for  transfer  of  vibrational  energy.  The 

ratio  7  =  <AE>a | | /<AE>d  has  been  evaluated.  A  reduced  energy  transfer 

SL  EXP 

quantity,  E^  =  <AE>d/RTe,  was  introduced  so  that  7  and  7  can  be 

calculated  from  universal  functions  for  E1  .  These  functions  exhibit  different 

e 

limiting  forms  for  the  EXP  and  SL  models. 

Comparison  with  exact  calculations  has  been  made  for  the  present  use  of  a 
fitted  classical  approximation  and  for  two  earlier  treatments  in  the 
literature  that  employed  the  WR  approximation.  Although  the  present 
approximation  is  relatively  more  accurate,  in  general,  the  treatment  of  Barker 
and  Golden  is  very  adequate  for  many  systems  and  conditions  and  is  simpler  to 
use.  Actually,  exact  calculations  are  not  arduous  and  are  also  advocated. 
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Peculiar  accidental  values 


Appendix 


It  seems  useful  to  summarize  the  nature  of  the  approximations  Involved  In 


the  T  and  BG  formulations.  Equation  7  with  E=EQ  can  be  written  as 

<  ^  V'.'  -v  (A1) 

x  v  '  -  '  T 

where  -  cl  =  -  ^ 

c2=  ^  AO  -  A v / 


C3=  2T  cr>- 

where  0(aE)=  J, p(  )  fy  (—■>  J 

The  sign  change  In  < aE> a i  -j  occurs  when  cl=c2*  l.e.  when  <^E>(j=<^E>1 


Troe's  expression  for  Fg  was  developed  to  take  Into  account  the  energy 
dependence  of  the  density  of  states  ratio  In  parameterizing  thermal  unlmolecular 
rate  constants  In  the  second  order  region.  As  a  result  there  Is  very  little 
resemblance  to  the  form  given  In  equation  (Al);  the  competition  between  "up" 
and  "down"  transitions  and  the  details  of  the  "down"  probabilities  are  absent. 


Using  the  Whltten-Rablnovltch  approximation  the  density  ratio,  D(aE)  can  be 
simplified  by  asumlng  a(E)  Is  a  constant  Independent  of  energy.  With  just  this 


constancy  assumption  D(aE)  Is  given  by  A2; 


u  iwm  t/  \  i  veil  wjr  n^|  i 

P„C&£)  -  I  +  1 


by  further  assuming  that 


-  A 


as  done  by  BG  equation  A3  results: 

n  .  .'*£  =  c*p  (('-•'}■ *=/ '(%  )  £))  (A3) 

H  ’V  ^ 

Expression  A 2  produces  Dq(aE)  <  D(A E)  for  all  AE;  the  difference  Increases  with 
Increasing  AE  and/or  molecular  complexity.  On  the  other  hand  when  A3  Is  used 


Dbg^  E)  <  0(  E)  for  small  AE  and  for  large  AE  Dgg  >  D(AE);  the  crossover  point 
Is  aE  ~  1700  cm"*  for  methyl  Isocyanide  and  Increases  with  molecular  complexity. 


Expression  A3  Is  a  better  approximation  to  the  true  density  ratio  for  all 

values  of  AE.  Thus  the  approximation  for  the  constancy  of  a(E)  which 

underestimates  D(aE)  Is  compensated  by  the  over  estimate  In  using  only  the 

linear  term  for  the  expansion  of  ln(l  +  x). 

The  approximations  used  for  D(AE)  will  determine  the  value  of  <AE>al  1  ( 

Al)  andy^XP;  Table  A1  Is  a  summary  of  calcul atlonal  results  for  representative 

EXP 

systems#  step  sizes  and  temperatures.  As  expected  from  the  form  of  FE»  y  j 
shows  the  largest  deviation  from  y^xp(exact).  The  Inadequacy  of  A2  for  D(^E)  Is 
shown  by  comparing  y  ^  and  Y  A3»  for  Increasing  step  size  and/or  molecular 
complexity  the  difference  between  these  quantities  Increases.  The  difference 
between  the  Integral  approximation  by  BG  (y  qq)  and  the  summations  used  In 
equation  Al  ('(  /&)  Is  negligible  and  for  practical  purposes  can  be  Ignored. 

The  goodness  of  the  BG  approximations  as  Illustrated  In  Fig  3  depends  on 
molecular  complexity#  step  size  and  temperature  In  a  complex  manner;  In  general 
the  BG  approximation  becomes  Inadequate  for  small  values  of  E0’.  For  a  given 
Ee'  the  difference  Increases  with  molecular  complexity. 


FXP 

Table  Al:  Calculated  Values  of  y 


<ae>.  t 

System  (cm—1)  1K1 

CH-.NC 

J  100  250 

500 
1000 
1500 
2000 

1600  250 

500 
1000 
1500 
2000 

CHT 

100  250 

500 
1000 
1500 
2000 

1600  250 

500 
1000 
1500 
2000 


EXPa  EXPb 

Y( exact)  !l _ 


-0.343 

-0.345 

-0.191 

-0.193 

-0.084 

-0.088 

-0.042 

-0.047 

-0.020 

-0.027 

-0.883 

-0.894 

-0.786 

-0.793 

-0.598 

-0.606 

-0.428 

-0.443 

-0.274 

-0.308 

-0.319 

-0.325 

-0.154 

-0.163 

-0.037 

-0.050 

0.009 

-0.012 

0.034 

-0.001 

-0.873 

-0.885 

-0.743 

-0.756 

-0.405 

-0.458 

0.037 

-0.161 

0.578 

-0.020 

EXPC 

EXPd 

BG 

YA3 

-0.345 

-0.344 

-0.192 

-0.192 

-0.086 

-0.086 

-0.044 

-0.044 

-0.022 

-0.022 

-0.894 

-0.883 

-0.792 

-0.787 

-0.601 

-0.599 

-0.424 

-0.423 

-0.261 

-0.261 

-0.325 

-0.324 

-0.162 

-0.161 

-0.047 

-0.047 

-0.001 

-0.001 

0.023 

0.023 

-0.885 

-0.875 

-0.755 

-0.751 

-0.440 

-0.439 

-0.017 

-0.017 

0.579 

0.578 

a)  Exact  y  using  calculated  from  equation  (Al). 

b)  y  calculated  using  equation  9. 

c)  y  calculated  using  equation  10. 

d)  y  calculated  using  equations  (Al)  and  (A3)  [constancy  of  a(E)  and 


expansion  for  £n  (1+x)]. 
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-0.344 

-0.192 

-0.086 

-0.044 

-0.022 

-0.883 

-0.783 

-0.603 

-0.436 

-0.287 


-0.324 

-0.162 

-0.047 

-0.001 

0.023 

-0.875 

-0.752 

-0.451 

-0.082 

0.353 


1  Inear 


e)  y  calculated  using  equations  (Al)  and  (A2)  [constancy  of  a(E)] 


Figure  Captions 


ci  cyp 

Fig  1:  Plots  of  7  (exact)  [solid  line]  and  7y^  (eq.  8a  [broken  line])  vs  E1 
TR 

(=<AE>^/RT  )  for  nitryl  chloride,  methyl  isocyanide,  cyclopropane,  and 

cycloheptatr iene  with  step  sizes  of  100,  200,  400,  800  and  1600  cm-*, 

each  at  250  K,  500  K,  1000  K,  2000  K  and  4000  K.  Equation  lib  in  text 

TR 

was  used  to  calculate  Tg  .  For  clarity  of  viewing,  only  data  points 

for  CH^NC  (circles)  and  eye  I oheptatr i ene  (triangles)  are  shown  for 

EXP  SL  SL 

7  (exact).  All  calculated  points  for  7-^  lie  on  the  7  (exact) 

curve. 

EXP  SL 

Fig  2:  Plots  of  7-j.p  /  7yp  [solid  line]  and  data  points  for 

FXP  qi  .  TR 

7  (exact  /  7  (exact)  vs  E^  (=  <AE>^/RT  )  for  reactants  and 

EXP  SL 

conditions  specified  for  Fig.  1;  where  7^  /  7^  = 

(E^/(l*E^)/((l-exp(-E^))/(l+exp(-E^)))  as  derived  from  eqs.  6  and  8a. 

Fig  3:  Plots  of  R^  [=  1(7^*  -  7^^(exact))  /  7^^(exact)l]  vs 

TR 

E^  (  =  <AE>^/RTe  )  for  CH^NC  (circles)  and  eye  I oheptatr iene  (triangles) 

with  step  sizes  of  100  (broken  line)  and  1600  (solid  line)  cm  *  using 

TR,  BG  and  T  approximations  for  the  EXP  model.  For  all  <  0.001  the 

data  points  are  placed  on  the  abscissa  (see  text).  The  filled  symbols  are 

for  R  >0  while  the  unfilled  symbols  are  for  R  <0.  The  difference 

Y  Y 

between  filled  and  unfilled  symbols  with  the  same  E'g  is  larger  than  what 

appears  on  the  plot.  Nonetheless,  the  absolute  errors  (comparison  with 

the  exact  calculations)  are  correctly  represented. 

EXP 

Fig  4:  Plots  of  7  (exact)  vs  temperature  for  methyl  isocyanide  (circles)  and 

cycloheptatriene  (triangles)  with  <aE>^  of  100  (filled  symbols)  and  1600 

/ _  .  ..  1  _  \  _  -  1 
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